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Accurate,  controlled  testing  of  high  technology  jet  pngines  requires  a 
fixed  installation  vhere  such  engines  nay  be  operated  throughout  their  full 
thrust  envelope  under  conditions  approxinating  the  installed  situation,  ar»d 
with  sufficient  instruaentation  to  assess  perfomance  paraneters.  These  in- 
stallations, jet  engine  test  cells,  have  taken  many  forms,  aaong  then  the 
"Hush-Souse"  installed  at  NAS  Hirasar,  CA  [Ref.  1]  fcr  installed  engine 
testing,  and  a Co an da  design  [Ref.  2]  for  noise  suppression.  The  most  common 
cell  design  is  typically  a block  house  type  installation  with  vertical  inlet 
and  exhaust  stacks  configured  for  velocity  profile  and  noise  control.  Arrange- 
ments with  horizontal  inlet  and  exhaust  stacks  exist  and  provide  more  uniform 
flow  profiles  to  the  engine  inlet,  but  the  large  dear  areas  needed  adjacent 
to  this  type  cell  precludes  its  frequent  use.  Mounting  hardware  and  monitoring 
equipment  are  provided  as  appropriate  to  the  engine  under  test. 

Pollution  control  is  a major  problem  in  the  operation  of  cells  with  today's 
high  power,  high  mass  flow  engines.  Noise  and  visible  and  invisible  pollutants 
are  emitted  in  quantity.  Judicial  action  initiated  by  the  State  of  California 
against  Naval  A-j.r  Station  facilities  in  California  has  brought  these  problems 
into  prominence.  Atmospheric  pollution  has  been  attacked  osing  various  forms 
of  water  droplet  adhesion,  mechanical  grid  entrapment,  or  electronic  ionization, 
with  baffle  combinations  for  acoustic  treatment  [Ref.  3].  Good  results  have 
been  obtained  at  NAS  Jacksonville,  FL,  with  a water  scrubbing  technique  to  re- 
move pollutant  particles.  All  these  methods  are  expensive  and  complicated, 
however,  and  simpler  ones  are  desired.  To  find  them  will  require  a detailed 
understanding  of  test  cell  aerodynamics. 
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In  a basic  arrangement  the  engine  is  situated  somewhere  near  the  center 
of  the  cell  to  alloy  a near  uniform  flow  profile  to  develop.  It  exhausts  into 
an  augmentor  tube  where  more  air  is  entrained  from  within  the  cell,  causing 
air  to  flow  around  the  engine  into  the  augmentor  tube.  The  ratio  of  entrained 
air  mass  flow  rate  to  engine  mass  flow  rate  is  known  as  "augmentation  ratio." 
The  entrained  "secondary  air"  acts  as  a coolant  for  the  hot  engine  gases,  ex- 
tending exhaust  stack  life,  and  as  a diluent  for  the  exhaust  products. 

Augmentor  parameters  such  as  engine-augnentor  spacing  and  augmentor  diameter 
are  important  to  proper  cell /engine  operation  due  to  their  effect  on  augmenta- 
tion ratio.  An  excessively  high  augmentation  ratio  (high  secondary  air  flow) 
may  cause  large  pressure  gradients  between  engine  inlet  and  exhaust  planes 
with  resulting  inaccurate  performance  measurements.  In  addition,  test  cell 
structural  limits  night  be  exceeded  due  to  excessive  cell  pressure  reduction. 
Insufficient  secondary  air  could  allow  hot  exhaust  gas  recirculation  to  the 
engine  inlet  with  performance  degradation  as  well  as  hot  spots  in  the  augmentor 
tube  and  exhaust  stacks.  Insufficient  secondary  air  flow  also  causes  excessive 
density  of  visible  emit  tan  ts.  This  say  violate  Ringlenan  Number  limitations  of 
local  pollution  ordinances  even  when  the  pollutant  output  quantities  are  within 
specified  limitations. 

Testing  for  optimization  of  engine-augnentor  relationships  cannot  easily 
be  accomplished  in  existing,  operational  test  cells.  Cells  are  scheduled  to 
narimm  capacity  for  economy  reasons,  making  test  runs  inconvenient.  Also, 
today's  large  pnginpg  and  soaring  fuel  costs  make  full  scale  testing  pro- 
hibitively expensive.  Clearly,  computer  modeling  of  flows  provides  a possible 
alternative. 

A working,  proven  model  would  be  able  to  predict  augmentation  ratios  and 
recirculation  of  hot  engine  exhaust  gases  *»s  a function  of  augmentor  design. 
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It  could  provide  information  on  the  optimal  locations  for  water  quenching 
devices  needed  to  cool  hot  afterburner  exhaust  flows  and  jet  break-up  devices 
used  for  noise  suppression  purposes.  These  capabilities  would  reduce  the  re- 
quirement for  expensive,  time  consuming,  trial  and  error  procedures  now  used. 

A number  of  models  with  some  or  all  of  these  capabilities  are  in  existence. 
One  by  United  Research  Laboratories  [Ref.  4]  uses  a one-dimensional  idealization 
combined  with  empirical  correlation  factors  to  adjust  ouput  values  of  exit 
static  and  total  pressure  (to  agree  with  experimental,  results),  then  generates 
performance  data  for  various  augmentor  combination.  £1 iin  and  Pucci  [Ref.  5] 
attacked  a similar  problem  in  modeling  a gas  eductor  system  for  gas  turbine 
powered  ships.  Due  to  differences  in  application,  however,  this  so  del  made  no 
allowances  for  spacing  between  nozzle  exhaust  and  tube  inlet,  a conaaonly  found 
configuration  in  test  cells.  Bailey  [Ref.  6]  also  used  the  one-dimensional 
idealization  and  accounted  for  friction  and  inlet  losses  through  empirical  data 
and  ingine-augmentor  spacing  by  means  of  theoretical  spreading  equations  [Ref,  7] 
for  the  engine  exhaust  stream.  Basinger  [Ref.  8]  also  utilized  a one-dimensional 
model  to  perform  a study  of  aircraft  ejector  optimization.  The  problem  was  also 
approached  using  two-dimensional  theory  by  Sayes  and  Netzer  [Ref.  9]  and  by 
Croft  and  Li 1 ley  [Ref.  10], 

Experimental  validation  of  the  1-D  model  of  Bailey  [Ref.  6J  and  the  2-D 
model  of  Hayes  and  Netzer  [Ref.  9]  was  the  subject  of  this  study.  The  one- 
dimensional Bailey  model  is  simple  and  uses  little  computer  time.  It  relies 
heavily  on  empirical  loss  factors  which  are  applied  to  flow  conditions  differ- 
ent from  the  original  experiments.  It  is  fundamentally  restricted  to  calcula- 
tion cf  trends  of  augmentation  ratio  with  variations  in  engine  flew  rate, 
engine  augmentor  spacing,  engine  diameter,  and  augmentor  diameter.  It  has, 
on  the  other  hand,  no  restrictions  to  conditions  of  application.  It 
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inccrpo rates  corrections  for  choked  flow  at  the  engine  nozzle  and  the  basic 
one-dinensional  equations  for  mass,  energy,  and  momentum  conservation  are 
not  limited  by  Mach  number.  It  does  not,  however,  incorporate  provisions  for 
handling  exhaust  shocks,  or  jet  spreading  for  choked  flow  conditions. 

The  two-dimensional  model  supplies  a far  more  complete  data  set.  It  pre- 
dicts temperatures,  velocities,  pressures,  and  turbulence  levels  throughout 
the  flow  field  within  the  augnentor  tube  and  parts  of  the  test  cell.  It  is  a 
more  rigorous  solution  of  the  basic  conservation  equations  and  relies  on 
empirical  constants  only  for  the  turbulence  modeling  and  boundary  conditions 
on  vorticity . It  is  limited  by  an  assumption  of  incompressible  flow  (more 
accurately,  density  is  assumed  to  vary  with  temperature  and  composition  but 
not  with  pressure)  which  prevents  accurate  results  above  flow  Mach  numbers  of 
about  0.6.  It  is  additionally  weakened  by  the  use  of  stream  function  and 
vorticity  as  primary  variables.  These  variables  make  the  solution  simpler  by 
eliminating  pressure  from  the  equations  and  by  reducing  the  dependence  on 
velocities.  However,  the  recovery  of  pressure  from  the  solution  is  extremely 
sensitive  to  calculations  for  the  stream  functions.  Both  these  areas  can  be 
strengthened  by  new  computation  techniques.  Tnese  include  elliptic  equations 
which  reduce  to  parabolic  in  appropriate  conditions  and  the  use  of  pressure 
and  temperature  as  primary  variables.  Tnese  techniques  are  currently  being 
investigated  at  the  Naval  Postgraduate  School. 

All  test  cell  models  in  existence  which  are  available  in  the  open  litera- 
ture have  the  same  basic  drawback  — none  have  been  experimentally  validated 
over  the  normal  range  of  test  conditions  found  in  actual  test  cells.  Due  to 
the  wide  range  of  assumptions  incorporated  in  the  models,  and  the  unusual  flow 
conditions  found  in  jet  engine  exhausts  (hot,  high  velocity  core  eni raining 
cold,  stationary  air),  models  must  be  tested  against  experimental  data  before 
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they  can  be  used  with  any  confidence.  Testing  should  include  measurement  of 
T£^oc^£y  profiles  in  the  test  cell  and  augaentor  tube  to  test  the  validity 
of  inlet  and  exit  profile  assumptions  cade  in  all  models  and  the  accuracy  of 
flow  velocity  computations  in  the  two— di-npncirmai  models.  Accurate 

measurements  of  flowrates  and  augmentation  ratio  are  necessary  to  provide 
inputs  to  the  2-D  model  and  to  check  predicted  results  in  the  1-D  model. 
Pressures  and  temperatures  along  the  augaentor  tube  and  in  the  test  cell  must 
be  measured  to  provide  data  for  comparison  to  2-D  predictions. 

To  make  all  these  measurements  in  a fall  scale  operational  r*>ri  would 
be  time  consuming  and  costly  for  the  same  reasons  mentioned  earlier.  For 
these  reasons,  the  subscale  test  cell  designed  and  built  at  the  Havai  Post- 
graduate School  [Ref.  11]  was  chosen  for  validation  efforts.  It  is  substan- 
tially less  expensive  and  more  convenient  to  operate  asd  can  easily  and 
inexpensively  be  completely  instrumented.  Configurations  ran  no  nn> /«r» j 
varied  and  data  easily  and  rapidly  collected  and  reduced. 


II.  KEIHOD  OF  PJVESTIGATIGN 


A previously  constructed  one-eighth  scale  turbojet  test  cell  [Ref.  11] 
vas  codified  to  increase  the  quantity  of  obtainable  data,  and  used  in  a 
valida'  icn  study  of  two  test  cell  flow  no  dels . A one-diaensional  and  a 


two-dimensional  model  were  evaluated.  The  cel  1 vas  used  to  determine  the 
effects  of  eagine-augnen tor  spacing,  engine  =ass  flow  rate,  engine  nozzle 
total  temperature  and  augaentor  inlet  geometry  on  augmentation  ratio  and 
augaentor  pressure,  velocity,  and  temperature  distributions . Experimentally 
measured  pressure,  temperature,  and  velocity  profiles  were  compared  with 
theoretical  predictions. 
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Ill . EXPERIMENTAL  APPARATUS 


A.  TEST  CELL 

Design  and  construction  of  the  sub  scale  turbojet  test  cell  are  detailed 
in  Ref.  11  and  shown  in  Figs.  1 and  2.  The  cell  is  a one-eighth  scale 
of  a HAS  Alary  da  test  cell.  A TF41  engine  was  scaled  to  one— eighth  in 
diameter  resulting  in  mass  flow  being  scaled  by  1/64  to  flow 

velocities  equal  to  those  in  the  full  scale  cell.  Air  was  drawn  into  the 
test  section  through  a horizontal  inlet  with  square  bellmouth  and  a flow 
straightening  section.  Tha  test  section  was  enclosed  by  hinged  pi  grig  less 
sides  to  allow  easy  access  and  vis  cal  monitoring  of  the  section  during 
operation.  The  augmentor  tube,  equipped  for  interchangeable  inlets,  exited 
the  cell  through  a removable  wall.  Its  downstream  end  was  attached  to  a 
deflector-plate-equipped  vertical  exhaust  stack.  The  stack  was  mounted  on  a 
wheel/rail  arrangement  which  allowed  translation  of  the  s tack/ augmen tor 
assembly  and  resulted  in  changes  in  spacing  between  the  engine  and  augmentor 
tube.  Removable  metal  grates  were  installed  in  tha  stack  to  permit  variation 
in  back  pressure. 

The  engine  used  to  simulate  turbojet/ turbofan  tailpipe  and  nozzle  condi- 
tions was  a forced  air  ramjet  supplied  with  c Depressed  air  free  an  Allis- 
Chalmers,  twelve-stage  axial  compressor.  Separate  three-inch  pressure  lines 
supplied  variable  combustor  (primary j and  bypass  (secondary)  air.  The  engine 
intake  was  simulated  by  a variable  suction  six-inch  line  drawing  through  a 
six-inch  bell  mouth.  The  combustor  was  a sudden  expansion  (or  dump)  type  fed 
by  nitrogen-pressurized  JP4  and  ignited  by  a methane-oxygen  torch  in  the 
combustor  wall.  The  arrangement  allowed  control  of  tailpipe  flow  rate, 
nozzle  total  temperature  and  pressure,  and  nozzle  geometry. 
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3.  ENGINE  AND  ADGMENIOR 

Figure  3 shows  the  basic  placement  of  the  ramjet  in  the  test  section  when 
viewed  from  downstream.  Figure  4 shows  the  geometry  of  the  engine-augaentor 
interface.  Hass  flow  rate  could  be  controlled  independently  in  the  primary 
and  secondary  lines,  providing  a simulation  of  a variable  bypass  turbofan,  if 
desired.  Mixing  of  the  primary  and  secondary  streams  was,  unfortunately,  not 
as  complete  as  in  an  actual  turbofan  and  in  hot  runs  the  arrangement  resulted 
in  a hot  inner  cere  of  primary  air  surrounded  and  pinched  by  cold  secondary- 
air.  Lengthening  of  the  mixing  section  will  eliminate  this  problem  in  the 
future,  Engi  ne-augaentor  spacing  was  varied  by  rolling  the  exhaust  stack 
along  its  rails.  The  three  augaentor  tube  inlets  primarily  tested  (Fig.  5) 
were  separate  pieces,  easily  changed  by  the  removal  of  two  screws. 

C.  INSTRUMENTATION 

Flow  rates  in  the  three-inch  primary  and  secondary  and  six-inch  suction 
lines  were  measured  using  ASHF.  type  flow  orifices  and  controlled  by  hand  valves. 
Pressure  taps  and  thermocouples  in  the  test  cell  allowed  monitoring  of  pressures 
at  the  engine  inlet  and  exhaust  planes  and  ambient  air  temperatures.  A pres- 
sure tap  and  iron-cons  tan  tan  thermocouple  in  the  combustor  cuter  wall  allowed 
calculation  of  nozzle  total  temperature  and  pressure.  Augaentor  flow  condi- 
tions were  monitored  by  27  pressure  taps  along  the  top  of  the  tube,  spaced  at 
one-inch  intervals  near  the  upstream  and  downstream  ends,  spreading  to  four 
inches  near  the  center.  Twelve  copper-constantan  thermocouples  were  also  placed 
along  the  tube  at  four-inch  intervals.  All  thermocouples  were  referenced  against 
an  ice  water  bath.  Velocity  profiles  in  the  augsentor  were  measured  with  a 
pitot  rake  consisting  of  seven  equally  spaced  small  diameter  total  pressure 
tubes.  The  rake  could  be  rotated  and  translated  to  obtain  vertical  and 
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horizontal  profiles  over  the  full  length  of  the  rube.  A thermocouple  was 
also  mounted  on  t!  s center  pitot  tube.  Velocity  profiles  in  the  cell  were 
measured  by  means  of  a direct  reading  anemometer  with  hot  wire  probe.  The 
probe  had  been  extended  sufficiently  to  allow  measurement  across  the  entire 
height  and  width  of  the  cell. 

D.  DATA  ACQUISITION 

The  data  acquisition  system  consisted  of  a fully  programmable  Hewlett- 
Packard  desk  top  calculator  with  a HP98673  Mass  Memory,  EP-9862A  Plotter,  and 
digital  tape  reader.  A3,  anc  F.  Model  ST 133  data  logger  scanned  45  pressure 


and  20  temperature  ch 


per  run  and  punched  raw  data  onto  a paper  tape 


for  reading  into  the  Hewlett-Packard  system.  The  data  reduction  system 
produced  reduced  flow  rates,  augmentation  ratio,  temperatures , and  pressures 
and  automatically  plotted  augmentor  pressure  and  temperature  profiles. 
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IV.  ElffERIMENTAL  PROCEDURES 

A.  TEST  PROCEDURES 

Initial  cell  startup  required  about  one  hour  for  warmup  of  the  Allis— 
Chalmers  compressor.  Once  air  was  directed  through  the  lines,  ten  to  fifteen 
minutes  were  required  for  temperature  stabilization  at  the  flow  orifices 
before  accurate  flow  rates  could  be  set.  Once  the  desired  mechanical  configu- 
ration of  the  cell  was  fixed,  flow  rates  were  adjusted.  Hough  adjustments 
were  made  using  water  manometers  located  at  the  hand  waives.  Find  adjustments 
were  made  by  running  flow  rate  data  through  the  data  reduction  system  and 
checking  calculated  values.  When  the  proper  values  were  attained,  a complete 
set  of  data  was  processed  by  setting  the  3.  & F.  data  logger  to  scan  the  40 
pressure  and  20  temperature  channels.  The  raw  data  was  recorded  on  punched 
tape.  The  tape  was  then  read  into  the  data  reduction  system  for  a complete 
data  readout.  When  cell  velocity  profiles  were  desired,  they  were  obtained  by 
manually  reading  and  recording  the  probe  measurements  as  it  was  inserted 
through  the  side  or  top  of  the  cell.  The  pitot  rake  was  also  repositioned 
manually  when  velocity  profiles  at  more  than  one  tube  position  were  measured. 

3.  MIA  RET-iETIOH 

Existing  programs  for  the  EP983QA  calculator  were  modified  to  haadla  in- 
creased data  quantity  and  extract  more  information  from  the  raw  data.  Flew 
races  were  calculated  using  ASHE  equations  for  &-.5D  orifice  pressure  tap  con- 
figurations. Temperatures  were  obtained  by  curve-fitting  published  thermocouple 
data.  Pitot  rake  pressures  were  combined  with  tube  wall  pressure  and  tempera- 
ture measurements  to  calculate  the  velocity  profile  at  each  position  of  the 
rake.  Inis  velocity  profile  was  then  integrated  by  Simpson’s  rule  to  find  the 
average  flow  velocity  in  the  tube.  The  velocity  profile  at  the  aft  end  of  the 
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augmentor  tube  was  quite  flat  and  varied  only  slightly  in  the  tangential 
direction.  For  all  tests  conducted,  the  average  velocity  at  this  location 
was  used  to  determine  the  flow  rate  through  the  augaentor  tube.  This  flow 
rate,  together  with  the  measured  engine  flow  rate,  allowed  calculation  of 
the  augmentation  ratio.  During  the  initial  phase  of  the  investigation,  the 
augmentation  ratio  determined  in  the  above  manner  was  validated  against  the 
value  determined  from  the  test  cell  velocity  profile.  A routine  was  included 
in  the  program  to  automatically  plot  tube  pressure  and  temperature  profiles 


when  desired 


C.  TEST  CONDITIONS 

Data  were  taken  for  the  set  of  test  conditions  presented  in  Table  1.  The 
bel  Inouth  data  were  used  for  1-D  Eodel  validation  and  the  flat  and  straight 
augsentor  inlet  data  were  used  for  2-D  node!  validation. 


TABT.F.  I 


TEST  CONDITIONS 


INLET 


Bellnouth 


Straight 


NOZZLE  FLOW  SATE 
(lb  /sec) 

Q 


ESG-AUG 

SPACING 


ENG-INLET 

SPACING 


-1.12D 


1.12D 


2.12D 


3.12D 


1.0 

1.0  (Hot) 

1.5 

1.0 

1.0  (Hot) 

1.5 

1.0 

1.0  (Hot) 

1.5 


B 


MODELS  TESTED 


A.  ONE-DIMENSIONAL  MODEL 

The  one-dinensional  so  del  tested  was  described  by  Bailey  in  Ref.  6. 
Geonetry  is  shown  in  Figure  6. 

The  model  assumptions  were: 


1. 

one-dinensional. 

steady 

-) 

X.  - 

adiabatic 

3. 

flow  uniform  at  cross  sections  1 

4. 

all  gases  ideal 

5. 

nixing  occurs  in 

a constant  area 

6. 

P » P 

1 2 

7. 

Aj^  + A2  * 

Values  of  tube  and  engine  nozzle  diameters,  nozzle  total  temperature  and 
pressure,  and  secondary  air  total  temperature  were  input.  A value  of  nozzle 
exit  static  pressure  (P^)  was  set,  which  set  ?2  , and  bypass  air  total  pressure 


(P  ) was  assumed  as  atmospheric.  Initially,  ?_  was  set  equal  to  ? 

x2  A?  * J2 
a “a 


Frog  these,  the  model  calculated  the  Mach  numbers,  and  velocities  at  1 and  2. 
The  loss  across  the  augsentor  inlet  was  then  calculated  using  an  empirical  loss 
constant  tor  the  particular  inlet  being  employed.  This  resulted  in  a new  value 


of  P and  the  process  was  repeated  until  successive  changes  in  P„  became 
z l2 


snail.  Primary  (1)  and  secondary  (2)  pass  flow  rates  were  calculated  and  used 
to  solve  the  one-dinensional  equations  of  nogen tua  and  energy  between  1/2  and  3. 

Spreading  of  a subsonic  jet  with  engine  augaentor  spacings  greater  than 
zero  was  handled  by  use  of  a spreading  parameter  to  determine  the  size  of  the 
nixing  zone  at  the  entrance  to  the  augaentor  tube.  Velocity  profiles  within 
the  nixing  zone  were  modeled  by  an  error  function.  The  ^-i-ring  zone  was  broken 
into  a large  number  of  snail  areas  and  then  properties  in  each  were  calculated 
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and  suszed  to  solve  the  one-dimensional  equations.  In  the  case  of  choked  flow 
(for  1-D  isentrooic  conditions)  in  the  nozzle,  a correction  for  choking  was 
applied  to  the  nozzle  flow  rate  and,  for  the  zero  spacing  case,  the  problem 
was  handled  as  before.  No  attempt  was  made  in  the  original  model  to  handle 
the  choked  flow  case  at  s pacings  greater  than  zero. 

Augmentor  tube  wall  friction  in  the  model  was  approximated  using  an 
equation  for  a flat  plate  drag  coefficient.  The  drag  coefficient  was  calcu- 
lated using  a Reynold * s Number  based  on  augmentor  length  and  V * (V^  + V^)/2 
and  T * (T.  + T^)/2  (for  determination  of  viscosity).  Although  not  mentioned 
in  the  text  of  Ref.  6,  an  empirical  coefficient  was  found  in  the  equation  in 
the  program  listing  which  arbitrarily  reduced  the  wall  friction  value.  The 
coefficient  proved  essential  to  successful  operation  of  the  model  since  without 
it,  the  calculated  flow  resistance  was  so  high  that  reverse  flow  was  predicted 
in  the  tube  for  all  cases  tested.  The  use  of  the  lower  value  of  wall  friction 
apparently  is  necessary  because  of  the  considerable  difference  in  flat  plate 
velocity  profiles  from  the  developing  coaxial  flow  within  the  augsentor  tube. 

To  use  the  model  as  written,  values  of  nozzle  total  temperature  and 
pressure  and  secondary  air  total  temperature  axe  input  with  geometries  and  the 
head  loss  factor  for  the  tube  inlet.  The  published  listing  was  modified 
slightly  to  allow  entering  the  experimentally  obtained  value  of  back  pressure 
(P^)  . In  the  original  model,  solutions  were  obtained  for  a number  of 
arbitrarily  set  values.  The  model  then  calculates  primary  and  secondary  flow 
rates  and  augmentation  ratio.  It  is  simple  to  use  and  requires  minimal  com- 
puting time.  At  the  same  time,  it  depends  heavily  on  empirical  constants  for 
inlet  losses  and  vail  friction  and  the  nplug  flow”  assumption  is  questionable 
in  the  case  of  high  velocity  entrance  flows.  It  provides  insight  only  into 
overall  behavior,  with  little  indication  of  behavior  within  the  tube  and  no 
information  on  flow  patterns  in  the  test  cell  itself. 
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The  no del  has  previously  been  tested  in  a l/24th  scale  plexiglass  cell 
[Refs.  6 and  12].  Bailey  found  the  nodel  reasonably  predicted  the  effects  of 
nozzle  total  pressure  on  augmentation  ratio  and  the  variation  of  secondary 
nass  flow  with  primary  mass  flow.  The  region  of  choked  nozzle  flow  vas  not 
investigated  since  no  provision  vas  aade  to  compute  the  spread  of  a supersonic 
exhaust  stream.  He  discovered  that  published  head  loss  factors  caused  exces- 
sive inlet  losses  and  predictions  were  most  accurate  with  a loss  factor  of 
essentially  zero.  No  conclusive  results  were  obtained  for  tests  of  engine- 
augaentor  spacing  effects  on  augmentation  ratio  due  to  widely  scattered 
experimental  data. 


B.  TWO-DIMENSIONAL  MODEL 


The  two-dimensional  nodel  tested  vas  described  by  Hayes  and  Netzer  in 
Ref.  9.  The  model  solved  the  basic  equations  of  mass,  momentum,  and  energy  in 
two-dimensional,  elliptic  fora.  Stream  function  and  vorticity  were  the  primary 
variables  for  mass  and  momentum  effects  which  reduced  coupling  problems  and  made 
equation  solution  easier.  Temperature  vas  the  dependent  variable  for  energy 
conservation.  Kinetic  heating  vas  ignored  and  specific  heat  vas  considered 
constant.  Turbulence  effects  were  included  in  the  form  of  a two  parameter 
model  describing  the  effective  viscosity.  The  relationship  used  vas: 

eit  y e 


where: 


y^ff  = effective  viscosity 
=“  empirical  coefficient 
p =■  local  density 
X = turbulence  kinetic  energy 
£ * turbulence  energy  dissipation  rate 
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The  equations  for  K and  e were  written  in  elliptic  form  and  the  five 
equations  were  solved  with  the  perfect  gas  law  to  describe  the  flow  field. 

Once  the  solution  converged,  axial  and  radial  pressure  distributions  were 
obtained  from  the  stream  function  distribution. 

The  equations  were  solved  using  point-by-point  Gauss-Seidel  with  tinder- 
relaxation.  The  model  was  assumed  axi symmetric,  and  described  by  a 43  by  40 
grid  in  cylindrical  coordinates.  A flat  augmentor  inlet  lip  (Fig.  5)  was  in- 
cluded in  the  model  as  an  approximation  to  the  bevelled  inlets  found  in  some 
cells  at  NARF  Alameda.  The  straight  pipe  was  approximated  by  raglting  this  lip 
very  snail. 

The  model  had  several  limitations.  The  elliptic  equations  using  stream 
function  and  vorticity  are  inherently  limited  to  low  subsonic  speeds.  In 
addition,  the  use  of  stream  function  and  vorticity  as  primary  variables  made 
calculated  pressures  extremely  sensitive  to  small  errors  in  convergence.  The 
rectangular  grid  system  does  not  allow  modeling  of  rounded,  bell  mouth  type  in- 
lets conoonly  found  in  operation.  Solution  required  an  average  of  170  minutes 
CPU  tine  on  an  IBM  360-67  computer.  On  the  plus  side,  it  is  a potentially 
valuable  tool  for  predicting  effects  of  engine-augmentor  spacing,  augmentor 
diameter,  and  cell  geometry  on  pressure,  velocity,  and  temperature  distribu- 
tions within  the  test  cell  and  augmentor  at  low  engine  flew  settings.  This 
information  is  needed  for  optimization  of  noise  suppression  and  chemical 
pollution  abatement  techniques. 
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VI.  RESULTS  AND  DISCUSSION 


A.  ONE-DIMENSIONAL  MODEL 

As  discussed  above,  the  no del  developed  by  Bailey  was  not  applicable  for 
choked  nozzle  flow  at  other  than  zero  engine-augmentor  spacing.  In  the  no  del, 
choking  was  assured  to  occur  for  the  one-dinens ional  isentropic  pressure  ratio 
(Y  * 1.4)  of  1.89.  The  flow  in  converging  no::zles  with  25°  half  angles  (as 
used  in  these  experiments)  does  not  actually  become  independent  of  back  pressure 
until  nozzle  pressure  ratios  in  excess  of  2.5  [Ref.  13]  are  obtained.  This  re- 
sults from  three-dimensional  effects  (radial  momentum,  etc.)  which  cause  the 
lines  of  constant  Mach  number  to  be  convex  in  shape  as  the  flow  exits  the  nozzle . 
In  the  experiments  conducted  in  this  investigation,  pressure  ratios  did  not  ex- 
ceed 2.4.  For  these  pressure  ratios,  the  flow  at  the  nozzle  exit  plane  would 
actually  be  slightly  subsonic.  The  flow  would  continue  to  accelerate  to  sonic 
or  slightly  supersonic  flow  conditions  with  resulting  weak  expansion  and  com- 
pression waves.  Since  the  expansion  and  compression  waves  are  quite  weak  for 
these  low  pressure  ratios,  a good  approximation  for  the  jet  behavior  may  be  to 
neglect  them  entirely  and  assume  that  the  jet  continues  to  accelerate  until  the 
jet  pressure  equals  the  local  ambient  pressure  a short  distance  from  the  nozzle 
exit  plane.  The  jet  behavior  could  then  be  approximated  by  considering  it  to 
spread  in  the  sane  manner  as  employed  by  Bailey  for  the  subsonic  jet.  This 
modification  was  incorporated  into  the  1-D  model. 

The  original  model  generally  predicted  trends  in  cell  performance  cor- 
rectly, although  with  same  severe  limitations.  Effects  of  engine-augmentor 
spacing  on  augmentation  ratio  are  s.  *wn  in  Figs.  7 and  8 for  the  test  cell  and 
model,  respectively.  The  predicted  augmentation  ratios  were  consistently  lower 
than  the  experimental  values.  Experimentally,  cell  augmentation  ratio  was 


found  to  decrease  slightly  and  then  level  off  with  increased  spacing  (Fig.  7). 
The  level  off  at  increased  spacing  indicates  that  jet  spreading  probably  became 
negligible  at  a point  about  two  nozzle  diameters  downstream.  The  decrease 
found  at  zero  spacing  probably  resulted  from  flow  interference  by  the  mounting 
flange  of  the  engine  nozzle  which,  at  these  separations,  was  very  near  the 
tube  inlet.  The  model,  which  predicts  continued  jet  spread  for  any  distance 
from  the  nozzle,  did  not  predict  the  leveling  off  found  experimentally  and 
could  not  predict  the  decrease  caused  by  flange  blockage.  Limiting  the  jet 
spreading  to  that  which  occurs  at  a distance  of  1.5  nozzle  diameters  did 
raise  the  predicted  augmentation  ratio  for  large  engine-augnentor  spacirgs 
(fig.  8). 

The  model  correctly  predicted  the  effect  of  changing  flow  rate  on  augmen- 
tation ratio.  Figure  9 shows  that  augmentation  ratio  decreased  steadily  with 
increasing  flew  rate  in  the  test  cell  and  Fig.  10  shows  the  same  trend  for  the 
model. 

As  mentioned  above,  the  model  predictions  for  augmentation  were  con- 
sistently less  than  the  measured  values.  This  is  more  evident  in  Figures  11 
through  14  for  different  engine-augmentor  spacings.  The  prediction  error  was 
greater  for  the  lower  nozzle  flow  rates.  Eliminating  augmentor  inlst  losses 
had  only  small  effects  on  the  predicted  augmentation  ratio  for  zero  engine- 
augmentor  spacing  (Fig.  11) c Eliminating  wall  friction  raised  the  predicted 
values  to  agree  with  experiment  at  low  flow  rates  but  to  overpredict  at  high 
flow  rates.  The  latter  was  probably  due  to  the  blockage  to  the  secondary  flow 
provided  by  the  expanding  jet  exhaust  in  the  experiments. 

Additional  comparisons  are  presented  in  Fig.  14  for  a 3.12D  engine- 
augsector  spacing.  Again  the  basic  model  predicted  the  correct  trends  but 
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underpredicted  the  magnitude.  The  model  predictions  are  seen  to  be  very  sensi- 
tive to  the  experimental  input  values  of  augmentcr  exhaust  pressure  (which 
depend  upon  exhaust  stack  resistance).  A change  of  0.1?  in  back  pressure 
changes  augmentation  ratio  by  approximately  32.  Apparently  this  is  the  most 
important  parameter  in  the  model  for  obtaining  the  correct  magnitudes  of 
augmentation  ratio.  Limiting  the  jet  spreading  ->  that  which  existed  at  a 
distance  of  1.5  jet  diameters  also  increased  the  predicted  values  as  did  re- 
ducing wall  friction. 

In  general,  results  agreed  with  those  found  fay  Bailey  for  subsonic  flow 
and  showed  the  same  trends  when  the  original  assumptions  were  applied  to  condi- 
tions which  would  be  choked  in  1-D  calculations.  The  model  predicted  trends 
accurately  but  requires  improved  wall  friction  and  jet  spreading  calculations. 


3.  TWO-DIHENS IONAL  MODEL 

Predictions  using  this  model  were  generally  very  accurate  when  ic  was 


operated  within  its  design  limitations.  Elliptic  type  equations  using  the 

stream  f uncti nn— vor  ticity  variables  are  generally  unreliable  for  nigh  subsonic 

velocities,  since  density  is  not  considered  to  vary  with  pressure  (which  has 

been  removed  from  the  equations) . Model  flow  rates  of  1.0  lb  /sec.  resulted 

* n 

in  a nozzle  exit  Mach  of  about  0.57  in  cold  flow.  The  higher  flow  conditions 
at  1.5  lb  /sec.  produced  an  exit  Mach  of  about  0.86.  Difficulty  was  experience- 

*5. 

in  running  comparisons  of  the  1.0  flow  rate  when  high  nozzle  total  temperatures 
(hot  flow)  were  eapJoyed.  The  model  requires  nozzle  exit  static  temperature 
as  an  input  parameter.  An  accurate  value  for  this  temperature  uas  impossible 
to  determine  from  the  experiments  since  the  exhaust  jet  consisted  of  a hot 


imw  core  at  greater  than  2090eR  surrounded  by  a cool  blanket  of  near  ambient 


secondary  air.  After  some  unsatisfactory  trial  and  error,  an  average  exit 
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tt^cr’r.-re  of  900°R  was  asstoed  for  all  hot  runs.  This  corresponded  to  an 
exit  Hach  of  about  0.8. 

Velocity  profile  predictions  for  cold  flow  in  the  augsentor  with  a flat 
plate  orifice  inlet  were  in  good  agreement  with  experiment  (Figs.  15  through  18). 
The  model  does  sees  to  consistently  predict  a slightly  slower  flattening  of 
the  profile  than  actually  occurred  since  it  consistently  overestimated  peak 
velocities  far  downstream  and  underpredicted  dose  to  the  nozzle.  The  least 
accurate  prediction  was  at  the  intake  lip  for  the  zero  spacing  case  (Fig.  15). 
This  is  most  likely  due  to  the  extreme  blockage  of  the  flow  at  the  tube  lip  in 
the  experimental  setup.  Since  the  calculated  velocity  sust  be  extracted  from 


st re  an  function  values,  any  error  in  calculation  of  the 


in  this 


rapidly  converging  region  would  be  somewhat  amplified  in  the  velocity  calcula- 
tion. As  anticipated,  prediction  error  was  slightly  greater  in  the  high  flow 
rate  case  (rig.  18).  This  was  not  surprising  since  "pinching"  of  the  nozzle 
flow  caused  a predicted  centerline  acceleration  to  a flow  Hach  number  of  1.12, 
about  2.5  nozzle  diameters  down  the  augsentor.  Surprisingly,  in  this  case,  the 
model  predicted  peak  -velocity  at  the  tube  inlet  to  within  32. 

Pressure  predictions  (Figs.  19  through  27)  were  in  relatively  good  agreement 
with  experiment,  although  not  so  impressive  as  the  velocity  results.  Pressure 
profiles  were  calculated  along  5 different  axial  grid  lines  of  the  tube,  ranging 
in  position  from  near  centerline  to  about  aid-radius.  The  outermost  lines  were 
just  inside  the  edge  of  the  lip  of  the  flat  intake  (Ref.  9).  The  line  J=17 
(the  second  outermost)  was  found  most  accurate  for  cases  with  the  lip  installed 
while  J*18,  (slightly  farther  out)  better  fit  those  cases  with  no  inlet  (rigs.  26 
and  27).  The  better  accuracy  is  found  to  occur  along  axial  lines  that  penetrate 
the  upstream  flew  in  the  least  disturbed  (more  axial  flow)  region.  In  addition. 


accuracy  is  generally  good  at  the  low  flow  rates  but  degrades  markedly  at  higher 
mass  flows.  These  results  were  expected  since  the  need  to  extract  pressure  field 
information  from  the  streamline  solution  through  two  successive  approximate 
derivative  calculations  renders  pressures  extremely  sensitive  to  snail  errors  in 
the  streamlines.  This  sensitivity  is  especially  evident  in  the  figures  where 
the  pressure  drop  in  the  entrance  section  of  the  augcentor  tube  is  consistently 
over  predicted,  the  error  being  aore  significant  for  the  high  flow  rates. 

The  model  proved  ineffective  in  predicting  hot  flow  results . An  example 
(Fig.  25)  shows  the  standard  result  — a grossly  excessive  pressure  drop 
accompanied  by  insufficient  rise  to  catch  downstream  valves.  Analysis  of 
model  predictions  showed  the  nozzle  flow  had  accelerated  from  a Mach  of  0.8 
at  the  nozzle  exit  to  a value  greater  than  Mach  1.2  along  the  centerline,  a 
flow  region  where  the  model  is  not  valid. 

Hesults  were  similarly  poor  for  hot  flow  temperature  predictions,  (rig.  28) 
where  data  was  non-ulmensionalized  using  the  maximum  temperature  of  the  run  to 
compensate  for  the  recognized  inaccuracy  of  input  flow  temperature.  Vhile  the 
model  consistently  predicted  gradual  miring  with  maximum  tube  wall  temperature 
at  the  downstream  end,  the  experimental  data  consistently  showed  peak  temperature 
at  about  one- third  the  length  of  the  tube  from  the  unstress  end.  The  discrep- 
ancy in  these  temperature  profiles  is  probably  attributable  to  both  the  high 
flow  Mach  numbers  and  the  poorly  defined  experimental  nozzle  exit  temperature 
and  velocity  distributions  as  discussed  above.  In  the  cold  flcv  case  (Fig.  29), 
agreement  was  excellent. 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 
A.  ONE-DIMENSIONAL  MODEL 

The  model  tested  has  several  advantages  over  the  core  elaborate  2-D  version. 
It  used  little  computer  time  and  was,  therefore,  very  significantly  less  expen- 
sive to  operate.  It  also  was  applicable  over  a wider  range  of  conditions,  not 
being  inherently  limited  to  low  subsonic  speeds.  The  model  was  not  capable  of 
accurate  quantitative  predictions  in  any  flow  condition,  but,  with  an  assumption 
of  no  physical  choking  at  the  nozzle  exit  and  only  weak  external  expansion  and 
compression  waves,  it  predicted  trends  accurately  up  to  flew  Mach  numbers  of 
about  1.2.  The  model  can  be  made  to  closely  fit  experimental  data  by  adjustments 
of  empirical  constants  (wall  friction,  jet  spreading)  and  then  be  used  to  consis- 
tently predict  conditions  for  a given  installation.  Predictive  accuracy  was 
very  sensitive  to  augmentor  tube  eshaus'.  pressure  which  in  turn  is  a function  of 
flow  rate  and  exhaust  stack  resistance. 

The  heavy  dependence  on  empirical  data  in  the  model  weakens  its  potential 
value  as  a general  design  tool  since  required  constants  may  change  from  installa- 
tion to  installation.  Loss  factors  for  unusual  configurations,  such  as  the 
bevelled  augmentor  inlets  found  at  NAHF  Alameda,  could  only  be  determined  by 
additional  experimental  measurements.  The  subscale  test  cell  utilized  in  this 
project  can  potentially  be  used  to  alleviate  this  shortcoming  by  providing  an 
inexpensive  and  convenient  test  device  to  establish  data  for  virtually  any  con- 
figuration. Validation  of  at  least  some  selected  conditions  against  full-scale 
test  cell  data  will  be  required. 
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B.  TWO-DIMESSIDNAL  MODEL 

This  model  is  inherently  limited  to  low  subsonic  flow  conditions,  a situa- 
tion which  renders  it  unusable  in  its  present  fora  as  a design  tool  for  high 
thrust  and  afterburning  conditions.  Its  quantitative  accuracy,  in  its  limited 
flow  region,  is  very  good  when  predicting  velocities  in  the  augaentor  tube,  but 
somewhat  weaker  when  predicting  pressures  due  to  the  need  to  extract  then  from 
streaa  function  information.  The  model  used  a large  amount  of  computing  time 
but  provided  an  enormous  amount  of  detailed  information  on  cell  performance  as 
a function  of  cell  design  and  engine  operating  conditions. 

The  greatest  value  of  the  model  probably  lies  in  its  function  as  a stepping 
stone  to  more  advanced  2— D models.  Current  work  is  being  directed  toward  a 
model  which  utilizes  velocities  and  pressure  as  the  primary  variables.  This 
model  should  improve  both  the  accuracy  and  utility  of  the  2-D  model  and  should 
provide  reasonable  solutions  for  the  choked  flow  operating  conditions. 
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